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Q«^ ■ Abstract 

- ' A lattice Boltzmann model for amphiphilic fluid dynamics is presented. It is a ternary model, 

jrt ' hi that it conserves mass separately for each chemical species present (water, oil, amphiphile), 

and it maintains an orientational degree of freedom for the amphiphilic species. Moreover, it 
models fluid interactions at the microscopic level by introducing self-consistent forces between 

' O [ the particles, rather than by positing a Landau free energy functional. This combination of char- 

acteristics fills an important need in the hierarchy of models currently available for amphiphilic 
fluid dynamics, enabling efficient computer simulation and furnishing new theoretical insight. 
Several computational results obtained from this model are presented and compared to existing 
lattice-gas model results. In particular, it is noted that lamellar structures, which are precluded 

^^ . by the Peierls instability in two-dimensional systems with kinetic fluctuations, are not observed 

^ ' in lattice-gas models, but are easily found in the corresponding lattice Boltzmann models. This 

points out a striking difference in the phenomenology accessible to each type of model. 

Keywords: Amphiphilic fluids, fluid dynamics, lattice Boltzmann, lattice gas, microemulsions, 
Peierls instability 
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1 Introduction 

Amphiphilic fluids consist of two immiscible phases, together with an amphiphile species that helps 
them to mix. A typical example is a mixture of oil and water with the addition of detergent. The 
molecules of detergent often have ionic "heads" that are hydrophylic, and hydrocarbon "tails" that 
are hydrophobic, and so there is a strong energetic preference for these to reside on the oil/water 
interface. Moreover, the presence of these orientable particles on the interface gives it a bending or 
curvature energy. These characteristics give amphiphilic fluids an extremely rich and complicated 
phenomenology, even in equilibrium. A wide variety of structures-including spherical and wormlike 
micelles, emulsion droplets, sponge phases and lamellae-can self-assemble on space and time scales 
that are very large in comparison with those of the component molecules, effectively ruling out 
molecular dynamics as a suitable methodology for studying such behaviour ^. The addition of 
fluid dynamics and nontrivial rheology to such complicated equilibrium behaviour thus creates a 
very difficult problem indeed, but one with relevance to a wide variety of industrial, chemical and 
biological applications. These include oil recovery, pollution remediation, and the dynamics of 
vesicles and biomembranes, to name but a few. 

Hydrodynamic lattice-gas automata (LGA) were first applied to reduced descriptions of am- 
phiphilic fluids - in which the amphiphile species is modelled only by its effects on the surface 
properties of two immiscible phases [2]. To capture much of the above-described phenomenology, 
however, requires a ternary vector model, in which a separate amphiphile species possesses both 
translational and orientational degrees of freedom. Such models have been used to study the equi- 
librium properties of amphiphilic fluids for more than a decade now |31. A few years ago, it was 
demonstrated that hydrodynamic LGA also provide an effective means of coupling the Hamilto- 
nian of such vector models to hydrodynamic flow with conserved momentum, thereby providing 
a self-consistent treatment of the hydrodynamics of ternary amphiphilic fluids [1]. This model, 
which can also be regarded as an extension of Rothman and Keller's immiscible-fluid algorithm [3] 
to include amphiphilic species, has been used to study the dynamics of phase separation [0], the 
shear-induced sponge-to-lamellar phase transition [7j, and many other phenomena of amphiphilic 
fluid dynamics. 

In the current work, we create a model of amphiphilic fluid dynamics that is based on a lattice 
Boltzmann (LB) model, rather than a LGA. Both models are kinetic in nature, and both have been 
implemented in both two and three dimensions [S], but the LB model is based on a single-particle 
distribution function, whereas the LGA is truly particulate. Phenomena that depend critically on 
kinetic fluctuations, nucleation and other particle discreteness effects are thus well described by the 
LGA model, but not by the LB model. On the other hand, the LB model is perfectly adequate for 
the treatment of bulk hydrodynamic phenomena, is simpler from an algorithmic point of view, and 
is more computationally efficient. Unlike the LGA model, which is far more complicated in three 
dimensions than in two owing to its use of lattice-specific many-body scattering processes, the LB 
model is easily implemented on lattices of any dimension. 

In addition to introducing the new LB model, an important purpose of the present work is to 
understand its limitations relative to the LGA model. Along these lines, we shall focus on the self- 
assembly of lamellar phases in two and three dimensions. According to a very general theorem of 
Peierls [HIE], in any infinite system with underlying kinetic fluctuations, thermodynamic equilibria 
that are periodic along one spatial dimension are necessarily unstable in two dimensions, but may 
be stable in three dimensions. Since the LB model lacks natural kinetic fluctuations, it may be 
expected to exhibit lamellae in two dimensions, and we present computer simulation evidence of 
this here. This situation stands in stark contrast to that for the corresponding LGA model, for 
which two dimensional lamellae have never been observed, presumably due to the Peierls instability. 



2 Lattice Boltzmann Models of Hydrodynamics 

LB models come in two varieties: So-called "top-down" models jllj assume an appropriate form 
for the free-energy, and assume that its approach to equilibrium is governed by a Ginzburg-Landau 
or Cahn-Hilliard equation; a LB collision operator is constructed that gives rise to the desired 
evolution equations in the hydrodynamic limit. "Bottom-up" models |12j . by contrast, add particle 
interactions at the kinetic level, and so the form (indeed, even the existence) of the free energy, and 
the time-dependence of the approach to equilibrium are emergent properties of the large scale dy- 
namics. Previous LB models of amphiphilic fluid dynamics have all been "top-down" in spirit, with 
a reduced description using one J3] or two ^^ order parameters, but no surfactant orientational 
degrees of freedom. The present model, by contrast is ternary (it contains three species, although 
one can also consider any single or pairwise combination of components) and adopts a "bottom- 
up" approach, making it easier to compare to existing LGA vector models. More importantly, 
unlike either the LGA model or existing LB models, it offers complete control of the viscosities 
and molecular weights of the various species present. Indeed, the algorithm allows for any num- 
ber of different species, including multiple types of amphiphiles, thus in principle enabling one to 
study inter alia the role of co-surfactants in these systems. This generality is evidently beneficial 
for studying a wider range of physico-chemical phenomena in amphiphilic fluids. Moreover, it has 
previously been argued that the current approach offers advantages of accuracy in the resulting 
hydrodynamic description |15j . 

A standard LB system involving multiple species is usually represented by a set of finite differ- 
ence equations ^1 EI , 

<(x + Ci,t + l)-<(x,t) = rjf i = 0,l,...,6, (1) 

where nf{x,t) is the single-particle distribution function, indicating the amount of species a (in- 
dicating, e.g., oil, water or amphiphile), having velocity Cj, at site x on a D-dimensional lattice 
of coordination number 6, at time step t. The collision operator 0,^ represents the change in the 
single-particle distribution function due to collisions. In the absence of chemical reactions, this 
operator must conserve the mass of each species, 

J2m^^-=0, (2) 

i 

where nia is the molecular weight for species a. Likewise, in the absence of external forces, it must 
also conserve the momentum of all the species, 

^^m.cA" = 0. (3) 

In the present work, we suppose that the system's microscopic degrees of freedom are attached 
to a heat bath at fixed temperature T, so that there is no corresponding requirement for energy 
conservation. 

For convenience, a "BGK" form JH] is often chosen for the collision operator, 

Qr = --«-nr), (4) 

where the equilibrium distribution function is constructed to satisfy certain desiderata, including 
the required conservation conditions jl9| I2f)j . Most LB models adopt a polynomial form for this 
equilibrium distribution function, such as 



nr''=9^n" 



l + (3oc^-u + -Piici ■ u)^ - /3oy + -Pii^r ' u)^ " /?o y (Ci • u) 



(5) 



where n'^ = J2i^iy ^^'^ where the degeneracy constants gi and the constant (3q are specified 
enth'ely by the choice of lattice. For example, /3o = 3 for the so-called D3Q19 lattice [201 1^ in 
three dimensions with 19 states. In Eq. (jSJ we have defined 



^='fr^ (6) 

where p'^ = rricjn'^ and p°'u°" = mf,-J2iCinf. Such a form for the equlibrium distribution ensures 
that the local mass and momentum values are collisionally invariant |22]. Indeed, one can directly 
verify that 

$:< = $:<, va 

i i 

and 

(7) 

^^m^Cj< = ^^m<^Cjnf, (8) 

<y i <y i 

where 

<^< + Of = <--«-nr) (9) 

is the post-collision distribution. The choice of n'^'^'^ then results in Navier-Stokes hydrodynamics, 
where each component a has viscosity v" = (tq- — 1/2) //Jq. 

3 Immiscible Lattice Boltzmann Model 

Immiscibility of two or more species is modelled by introducing a discrete realization of the self- 
consistently generated mean-field body force between the particles |12j . Without an amphiphilic 
species, and considering only nearest-neighbour interactions among the immiscible species a G 
{1, . . . , am}, the simplest possible form for such a force is 

F'^'^(x,t) = -r(x,t)^^5,^V'^(x + c,,t)c,, (10) 

CT i 

where the parameter ^o-o- (positive for repulsive forces, negative for attractive ones) represents 
the interaction strength between the two components, a and a. Here we have introduced the 
effective charge ip'^ , which can be chosen to be any suitable function of the local mass density p"^, 
as appropriate for the force laws of the hydrodynamic system of interest. The dynamical effect of 
the force is realised in the BGK collision operator in Eq. (0J) by adding to the velocity u, in the 
equilibrium distribution nf "^ of Eq. Q, an increment 

6u'' = t^F'^'Vp^ (11) 

This procedure has recently been shown to be systematically derivable from the body-force term of 
the standard continuum Boltzmann equation, based on a Hermite moment expansion j^B]- Indeed, 
with this algorithm, phase separation has been demonstrated for sufficiently large values of g^a |121 
I24| . It is worthwhile to point out that, although the overall momentum change at any given lattice 
site and time may be non-zero, 

F(x,t) =5]F'^''=(x,t) =5]p'^(x,t)(5u'^(x,i), (12) 



Newton's third law, and hence global momentum conservation, is satisfied exactly by Eq. ()1U() . 

EF(x,t) = 0. (13) 

X 

4 Amphiphilic Lattice Boltzmann Model 

With the addition of amphiphile, the LB Eqs. (^, @, © and ()11|) remain unchanged, and the 
species label a simply extends to include the amphiphile species s. Although as noted above the ba- 
sic algorithm is general enough to treat multiple amphiphile species, for simplicity here we consider 
only a single type of amphiphile. Incorporation of the physics of amphiphilic fluid particles requires 
two fundamentally new types of body forces, in addition to those described above. Unlike ordi- 
nary species such as oil and water, which can be approximated by point (or spherical) particles so 
that their interactions depend on relative distances alone, the essential character of an amphiphilic 
molecule is more appropriately represented by a dipole; their interactions depend not only on their 
relative distances but also on their dipolar orientations |1]. It is these additional fundamental 
features which result in all the fascinating phenomenology described in the Introduction. 

We first introduce an average dipole vector d(x, t) at each site x to represent the orientation of 
any amphiphile present there. The direction of this dipole vector is allowed to vary continuously in 
the present model (although we could do so, we do not attempt to discretise the dipolar orientation 
here). We assume that the dipole vector for a single amphiphile molecule has magnitude do, and we 
suppose that there are A^ amphiphile molecules represented by a single lattice-gas particle, so that 
|d(x, t)| < do = Ndf), where equality holds if and only if all the amphiphile molecules are perfectly 
aligned. Note that the model does not specify dipole information for each velocity Cj, for reasons 
of computational efficiency and simplicity; information on the average dipole direction at each site 
is adequate for the present model. The propagation of d(x, t) can then be defined by the simple 
relation, 

n"(x,t + l)d(x,t+l) =^n|(x-Ci,t)d(x-Ci,t). (14) 

i 

Note that the dynamics of the flow is entirely contained in the distribution function hf (x — Cj , t) . 
In Eq. dJ, d represents the post-collision average dipole vector per site, whose relaxation is also 
governed by a BGK process, 

d(x, t) = d(x, t)- — [d(x, t) - d'=^(x, t)] (15) 

with T^ representing relaxation of a dipole vector to a local equilibrium d^'^{x,t). Eqs. H14|) and 
(|15)) completely determine the evolution of the dipoles once an explicit expression for d'''^(x, t) is 
specified. Furthermore, it can be seen that the magnitude of the dipole vector is guaranteed to be 
less than do at all times if r^^ > 1 and |d'^'5| < do- 

Unlike momentum vectors, there is no conservation requirement associated for the dipole vec- 
tors. On the other hand, their direction and magnitude are affected by the distribution of the other 
dipoles as well as of oil and water molecules in their neighbourhood. It is physically reasonable to 
assume that the local equilibrium at a given site is the same for all amphiphile molecules, regardless 
of their velocity values. That is, if we let d^^(x, t) be the equilibrium dipole value for amphiphiles 
of velocity Cj , then 

d^'^(x,t) = d^^(x,t), Vi. (16) 

Moreover, such an equilibrium property should not depend on the history of the dipole distribution. 
Both the orientation and magnitude of the equilibrium dipole vector are dictated by the nature 



of the particles in its neighbourhood and so one of most appropriate choices is to use the Gibbs 
measure: 



d°'i(x,t) = do- 



J2 e-l3Ha{x,t) 



(17) 



where 0„, for a S {1, ... , W} are a set of unit vectors representing all possible dipole orientations. 
Obviously, Eq. H17|) ensures |d®^| < do. The parameter /? = 1/T represents the inverse of the 
temperature of the aforementioned heat bath, with which the microscopic relaxation process is 
assumed to be in contact. From Eq. H17() . one can immediately make the interpretation that the 
dipole tends to align itself in the direction which minimizes the energy. Ha- If we consider the 
dipoles' vector interactions, it is plausible to assume the following form for Ha, which can be 
interpreted as a generalization of the Heisenberg model Hamiltonian in magnetism, 



Ha{:>c,t) = -4>^ •h(x,t) 



(18) 



where 



h(x,t) = h"(x,t) + h=(x,t) (19) 

is a mean field created by the surrounding distributions. Its two contributions, h^(x, t) and h'^(x, t), 
are associated with distributions of ordinary species and amphiphiles, respectively. If we consider 
only nearest-neighbour interactions for simplicity, these vector fields have the following forms 



h'^(x, t) =^e„^ 7i'^(x + Ci,t)ti 



(20) 



where e^ (which may, for example, take its values from the set {— e, 0, +e} where e is constant) is 
the colour charge (order parameter) for ordinary species. This form may be interpreted physically 
as a discrete approximation of the colour gradient of the immiscible species. Similarly, we have 



h^(x,t)=E 



^n|(x + Cj,i)0j •dj(x + Cj,t) + n- (x, t)dj(x, i) 



where we have defined the traceless second-rank tensor 



e, 



D^. 



(21) 



(22) 



and where in turn I is the second-rank unit tensor and D the spatial dimension. The physical 
reason for the appearance of the tensor Oj is that dipoles tend to align when their relative positions 
are orthogonal to their mean dipolar orientation (i.e., d • Cj = 0), and anti-align when these are 
parallel. Since we let 0„ take on continuous values, in three dimensions Eq. (|17() can be integrated 
analytically to give 

d"'^ = do coth iPh) ^ 



(5h, 



which becomes 



" 3 



(23) 



(24) 



in the limit of small (3h. Here we have defined the unit vector h = h(x, t)/|h(x, t)|, and magnitude 
h= |h(x,t)|. 

With the presence of an amphiphilic species s, the force on an oil or water particle a needs to 
include an additional term 

F'^(x, t) = F'^'^(x, t) + F'^'^(x, t), (25) 
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where F""''^ is just the original contribution from the ordinary molecules given by Eq. H1U|) : the 
second term in Eq. (|25|) is the contribution from nearby amphiphiles whose specific form can be 
derived directly from a leading-order expansion of Eq. (|10() . To see how this latter expression may 
be derived, recall that an amphiphilic molecule possesses a hydrophilic (positive colour) head and 
a hydrophobic (negative colour) tail. Such a colour dipole can be modelled by a pair of oppositely 
(colour) charged fictitious particles with locations displaced by ibd/2 from the molecular center 
of mass location x. On this basis, we can derive an analytical form for the resulting force from 
Eq. (|10j) via a Taylor expansion in d. When only nearest-neighbour interactions are considered, 
the form of this interaction is given by 






(26) 



where the parameter g^g is the coupling coefficient between an ordinary and an amphiphilic species. 
Using similar physical arguments, we can derive the forces on amphiphilic particles. 



F^(x,t) = F^'^(x,t) + F^'^(x,i), 



(27) 



where F"''^ and F'^''* represent forces from ordinary species and amphiphiles, respectively. In similar 
fashion, explicit forms for both of these forces can be derived directly from the force law, Eq. (|1U() . 
for ordinary immiscible species; physically, the force F*''^ is simply the reaction force to F*^'*. It 
may be shown to have the following form. 



F^''=(X, t) = 2V'^(X, t)d(x, t) ■Y.dasY. ^iV'"(x + Ci, t), 

from which we can verify that Newton's third law is satisfied, 



(28) 



E 



^F'^'^(x,t) + F^''=(x,t) 



0. 



(29) 



To the leading order of a Taylor expansion in the ratio of jcj| to the colour gradient scale length, 
Eqs. (Uni) and ((2HI) reduce to 



26c2 



D{D + 2) 



9.sr(x){v2[d(x)V^^(x 



DVV 



d(x)V^(x)]} 



and 



F"'^(x) 



26c2 



D{D + 2 



-5<xsV''(x)d(x) • ^ [lVV"(x) - DVVri^) 



(30) 



(31) 



The above expressions make manifest the physically reasonable result that fiuid particles experience 
no net force if the other particles in their neighorhood are distributed uniformly. More interesting 
is the result that the leading-order contributions involve second-order spatial derivatives. This 
implies that amphiphiles tend to reside near the interface between two immiscible fluids. 

Finally, F**'* represents the force due to interactions between amphiphilic particles. This de- 
pends not only on the relative distance but also on the relative orientation of the dipoles in question. 
Though considerably more involved algebrically, its expression is still derivable from small expan- 
sions of the polarization d about x as well as those at the neighbouring sites, as described above. 



In particular, when only nearest-neighbour interactions are considered, this force has the simple 
form 



+ 



d(x + Ci)d(x) + d(x)d(x + Ci)] • Ci} tfj'i^ + Ci). (32) 



Here the coupling coefficient gss should be negative if we wish to model attraction between two 
amphiphile heads and repulsion between a head and a tail. Likewise, it can be verified directly that 
global momentum conservation is satisfied by this form; that is, 

^F^'^(x,t)=0, (33) 

X 

when boundary effects are absent. 

To summarize, the basic LB evolution equation is given by Eq. (^. The BGK operator for 
collisional relaxation is given by Eq. @, where the collisional equilibrium is given by Eq. (jSJ; 
this describes an ordinary LB algorithm for miscible Navier-Stokes fluids. To realise immiscibility, 
we implement the force, Eq. ()1U() . between neighbouring particles by shifting the velocity used 
in Eq. Q by the amount specified in Eq. pijl : this describes the previously known "bottom-up" 
immiscible fiuid LB algorithm. The principal contribution of this work is the introduction of a 
surfactant species with an orientation that propagates according to Eq. (|14p . and relaxes according 
to the BGK-like orientational collision operator, Eq. H15|) . with orientational equilibrium given by 
Eqs. (|23|) . (|2fl|l and (|2T|) . The force on ordinary oil/water species due to amphiphile is given by 
Eq. ()26() . that on amphiphile due to ordinary species is given by Eq. H28() . and that on amphiphile 
due to other amphiphile is given by Eq. (|32|) . The effect of these forces on the collisional relaxation 
is also treated by shifting the velocity used in the BGK equilibrium by an amount specified in 
Eq. (jmi . Taken together, these equations completely specify the new ternary amphiphilic LB 
model and associated algorithm. 

5 Simulations 

In this section we demonstrate qualitatively that some of the basic phenomenology of amphiphilic 
fluids is indeed exhibited by the new model. We begin by verifying that the model with the 
amphiphilic interactions extinguished displays complete oil-water phase separation. All simulations 
shown in this section were performed on a domain of size 32 x 32 lattice units with periodic boundary 
conditions, and used To- = r^ = r^ = 1 and (3 = 10. 

Fig. n shows the equilibrated colour order parameter at 6000 and 16000 time steps after an off- 
critical quench. The run began with blue density n^ = 0.45, red density n^ = 0.65 and surfactant 
density n** = 0.25, distributed homogeneously in space, aside from a small random fluctuation. For 
sufficiently strong interaction parameter (g^r = 0.03 was used here) between the immiscible red and 
blue phases, complete separation was observed as can be seen in the figure. Because the amphiphile 
interaction parameters were turned off (^f^^ = gss = 0), the separation did not arrest, and the blue 
fluid formed a single droplet (recall that the boundary conditions are periodic). 

Next, using the same density values and initial conditions as above, but with the amphiphilic 
interactions switched on (gijs = gss = —0.01), Fig. |21 shows the colour order parameter and the 
surfactant director density after 16000 time steps, by which time the system appears to be in 
a stationary state. Note that the amphiphile has organised itself on the surface, with directors 
pointing from blue to red on average. This effectively demonstrates that the amphiphile has arrested 
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Figure 1: Two-dimensional binary oil(red)-water(blue) phase separation. Colour (oil- water) distri- 
butions shown at lattice timesteps t = 6000 (left) and t = 16000 (right) 

the phase separation; if the separation were to continue, there would not be sufficient surface area to 
accommodate the amphiphile, which incurs a large energy penalty for being away from a surface. 
Such arresting of domain growth has previously been observed in LGA studies of amphiphilic 
fluids m. 

Fig. ini shows the equilibrated colour order parameter and director density, 20000 time steps after 
a critical quench on a domain of size 32 x 32 lattice units. The run began with n'' = n^ = 0.5 and 
n^ = 0.1, distributed homogeneously in space, aside from a small random fluctuation. No further 
evolution of the order parameter was observed after a few thousand time steps, strengthening the 
case that this is indeed an equilibrium state. 



6 Discussion and Conclusions 

As mentioned in the Introduction, an important focus of our studies to date with the new model 
has centered on the question of the existence of lamellae in two spatial dimensions. The theorem 
that we invoke was derived by Peierls in the context of electronic structure |H1 El j but is general 
enough to apply here. Let us suppose that a structure (p{x) that is periodic along dimension x 
minimizes some free energy functional F[(p]. Let us further suppose that this structure undergoes 
thermally induced fluctuations in a Gibbsian distribution exp(—6F/T), where 5F is the free energy 
excess due to the fluctuation. Peierls demonstrated that the variance of the component of such 
fluctuations in the direction of the periodicity is proportional to 



dkx k±dk± 



(34) 
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Figure 2: Two-dimensional domain growth in a ternary (oil-water-surfactant) fluid. Colour (oil- 
water, left) and amphiphile director (right) distributions are both shown at lattice timestep t = 
16000 
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Figure 3: Formation of a lamellar phase in a ternary oil-water surfactant lattice Boltzmann fluid. 
Colour (left) and amphiphile (right) distribution are shown, both at lattice timestep t = 20000 
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where (p is an analytic function of its arguments that vanishes when its arguments are zero. This 
variance is easily seen to diverge logarithmically for small wavevector (large system size) for any 
finite temperature T > 0. Moreover, similar arguments show that such a divergence does not occur 
in three spatial dimensions. 

Thus, we should expect lamellar structures to be unstable in an infinite two dimensional system 
with finite-temperature fluctuations. Indeed, previous lattice-gas studies of amphiphilic fluids in 
two dimensions |3 [7j have not found stable lamellae, despite extensive parametric searches for 
them, except when the system was sheared jjj. Lattice-gas studies in three dimensions, however, 
do exhibit stable unsheared lamellar structures jH|- We suggest here that the observed instability 
of lamellae in two dimensions is due to the Peierls mechanism. However, we must qualify this 
claim with the proviso that it is somewhat dangerous to assume that the Peierls instability applies 
in a lattice-gas simulation: although a LGA is indeed a finite-temperature, fluctuating system, 
computer simulations of it always have a finite domain size, and hence a minimum sustainable 
wavevector. This objection might be mitigated if it could be shown that lamellar phases reappear 
at zero temperature. Unfortunately, it is impossible to lower the effective temperature of a LGA 
to zero; kinetic fluctuations will always be present by the very nature of the model. LB methods, 
by contrast, have no kinetic fluctuations. Thus, if an LB model were to exhibit lamellae in two 
dimensions, it would strengthen the argument that their absence in the LGA model is due to the 
Peierls instability. As described in the previous section, we tested this by simulating the critical 
quench shown in Fig. |31 that this resulted in a lamellar phase supports the above interpretation of 
the Peierls instability. 

In conclusion, we have described a ternary lattice Boltzmann model of amphiphilic fluids. This 
model differs from previous LB models in that it accounts for amphiphilic director orientation, and 
incorporates microscopic interactions between the fluid particles, without the need to postulate a 
free energy functional. It differs from LGA models in that it does not include fluctuations associated 
with particle discreteness. The lack of fluctuations gives the model a different phase behaviour from 
the corresponding LGA; in particular, we have presented evidence from computer simulation that 
it exhibits lamellar phases in two dimensions, which are prohibited by the Peierls instability in 
models with fluctuations. 
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